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ABSTRACT 

The  equations  of  motion  of  compressible  viscous  and  heat-conductive 
fluids  are  investigated  for  initial  boundary  value  problems  on  the  half  space 
and  on  the  exterior  domain  of  any  bounded  region.  The  global  solution  in  time 
is  proved  to  exist  uniquely  and  approach  the  stationary  state  as  t  +  **, 
provided  the  prescribed  initial  data  and  the  external  force  are  sufficiently 
small.  _ 


AMS  (MOS)  Subject  Classifications:  35Q99,  76N10 

Key  Words:  Equation  of  general  fluids ,  Initial  boundary  value  problem. 
Global  solution  in  time 

Work  Unit  Number  1  (Applied  Analysis) 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 . 


SIGNIFICANCE  AND  EXPLANATION 


The  smooth  motion  of  compressible  viscous  and  heat- conductive  fluids  is 
described  by  a  nonlinear  system  of  five  partial  differential  equations.  The 
first  mathematical  question  for  such  a  system  is  the  existence  and  uniqueness 
of  the  solution  for  the  initial  and  boundary  value  problem.  Since  the  local 
existence  in  time  of  the  unique  solution  is  known,  the  main  problem  is  the 
existence  of  the  global  solution. 

In  this  paper  we  show  the  existence  and  uniqueness  of  a  global  solution 
in  time  and  its  decay  to  the  stationary  state  as  time  tends  to  infinity  for 
the  particular  case  of  the  initial- boundary  value  problem  on  the  half  space 
and  on  the  exterior  domain  of  any  bounded  region,  provided  the  initial  data 
and  the  external  force  are  small.  The  analysis  which  uses  the  energy  method 
extends  our  previous  works  concerning  the  initial  (boundary)  value  problems  on 
the  whole  space  and  on  the  interior  domain  (MRC  Report  #2194,  2237);  some  new 
a  priori  estimates  for  the  domains  considered  here  are  required. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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INITIAL  BOUNDARY  VALUE  PROBLEMS  FOR  THE  EQUATIONS  OF  MOTION 
OF  COMPRESSIBLE  VISCOUS  AND  HEAT-CONDUCTIVE  FLUIDS 


Akitaka  Matsumura  and  Takaaki  Niahida 

|1.  INTRODUCTION 

The  action  of  viscous  compressible  fluids  is  dascribad  by  tha  system  of  five  aquations 
for  tha  danaity  p,  tha  valocity  u  «  (u1 ,u2.u3)  and  tha  taaparatura  6s 


Pt  ♦  <PuJ> 


(1.1)  u*  +  U3U*  +  £  p  -  ”  +  U3  )  +  u’v  v  fi* 

t  xj  p  xi  p  xj  xi  \  xi 


i  -  1,2,3, 


e  +  u3e  ♦  ~  u3  -  ~  {cce  >  ♦  »>  , 

t  xJ  PCV  pc7  Xj 


where  p  la  tha  preaaura,  y  is  tha  viscosity  coafficiant,  u'  is  the  second  coefficient 
of  viscosity,  <e  is  tha  coefficient  of  heat  conduction,  c?  is  the  specific  heat  at 
constant  volusw,  all  of  which  are  known  functions  of  P  and  0,  and  1  is  the 
dissipation  function! 

T  -  |  (u3  ♦  u*  >2  ♦  M'(u*  )2  . 

2  *k  j  Xk 

He  consider  the  initial  boundary  value  problem  of  (1.1)  in  the  region  t  >  0,  x  e  ft. 
The  boundary  condition  is  supposed 


(1.2) 

or 

(1.3) 


Ill  •  ul  -  0, 

8l 

- 

el 

1 3ft  1- 

1 

30 

1- 

38 

i!, 

!|  Ul  0, 

ho  l- 

3n 

|ao 

3n| 

t  >  0  , 


t  >  0  . 


The  initial  condition  is  given  by 

(1.4)  (p,u,8)(0,x)  -  ( Pq ,uq , 8g ) (*) ,  x  «  ft  . 
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The  local  existence  theorem  for  the  initial  boundary  value  problem  (1.1)- (1.4)  is  proved  by 
Tani  (10]  under  full  generality.  Here  we  want  to  solve  the  problem  (1.1)— (1.4)  globally  in 
time  under  the  following  assumptions. 

1,1.  The  domain  Q  is  the  half  space  R^  *  {x  e  R3,  x^  >  0}  or  an  exterior  domain  of  any 
bounded  region  with  smooth  boundary.  Here  we  note  that  the  initial  (-boundary)  value 
problems  on  the  whole  space  and  on  the  interior  domain  are  solved  globally  in  time  in  [6] , 
[7]  and  [6]. 

A. 2.  w, p' , x,p  and  are  smooth  functions  of  P,  8  >  0,  and  P,  K.c^p.p^.pg  >  0, 

u‘  +  |  V  >  0. 

4 

A. 3s  The  external  force  f  is  given  by  the  potential  #(x)  e  H  (ft), 

d.5)  f1  -  ,  i  -  1,2,3,  x  e  a  , 

xi 

where  and  in  what  follows  we  use  the  Sobolev  space 

a  a  a 

h*(B)  -(f«  l2(8),  Dkf  -  oVax^Sxj  8X3  ,  | a|  -  k}  e  l2(b>,  i  <  k  <  »} 

with  the  norm  If  I  »  (  J  /  |Dkf  I  2dx)1/^2. 

k-0  0 

A. 4.  The  initial  data  are  smooth  functions  close  to  a  constant  state  (p,0,"6),  where 
p  and  S  (also  appeared  in  (1.2))  are  any  positive  constantns,  i.e., 

PQ  -  P,  uQ ,  8q  -  0  e  H3  ( 8)  and 

(1.6) 

lpfl  -  p,  uQ,  9Q  -  "S »3  is  small  . 

A. 5.  The  compatibility  condition  on  the  initial  and  boundary  data  is  satisfied  as  follows  1 


u  “0,  9  «9 

0  an  0  an 


*  T0  'Vo*.,,  *  *  •<<«/’>., U  '  M; 


where  we  used  the  notation  *  p ( pQ , )  etc.  In  the  case  of  Neumann  condition  (1.3)  for 
6  we  need  necessary  change  of  compatibility  condition. 

He  use  the  following  notations  for  the  function  spaces: 

0^(5):  Banach  space  of  bounded  continuous  function  on  R, 

cSfl)  *  (fe  C°<5),  D^f  e  C°(fl),  1  <  k  <  i>, 

C^(t^ /t^/B)  ■  ju(t)  :  i-times  continuously  differentiable  function  of 

t  e  tt^.t  ]  with  values  in  a  Banach  space  B,  where  the  norm  is  given  by 

Ifo 

L^lt^t^tB)  “  (u(t)  :  L2~f unction  of  t  e  (t^ .t^l  with  values  in  B,  where  the 


max  sup  •)  I  }, 

0<k  <1  ^<t<t2 


•2-)1/2) 

The  stationary  solution  (p,u,8)(x)  of  in  a  neighborhood  of  (Pf0#*b)  in 


norm  is  given  by  (/  lu(t# •) l*dt }. 

s 


2  3 

H  (R)  has  the  form  (cf.  lemma  2.1) 


(1.8) 


p(x)  p  (n, 0)  _  ^ 

j  —~r —  dn  +  #(x)  »  o,  u  •  o,  9  -  8 


He  can  state  the  main 

Theorem  1,1,  under  the  assumptions  A.1  ~  A. 5  there  exiBtB  a  constant  e  >  0  such  that 
if  »P0  -  P.  V  60  “  **3  +  "'4  «  E  ' 

then  the  initial  boundary  value  problem  (1.1)-(1.4)  has  a  unique  solution  ( p,u,  8) 
globally  in  time  and  a  unique  stationary  state  ( P,0, 8)  which  satisfy 


(1.9) 

and 

(1.10) 


p  -  p  e  C0(0,-jh3(R)>  n  C1(o,-,h2(R)) 
u,e  -  5  e  C°(o,»jh3(0))  n  C1  (0,-ni1  (Q)) 

|p  -  p,u,8  -  "S(t)  I  n  ♦  0  as  t  ♦  •  . 

|c°(fl) 


The  theorem  is  proved  by  the  energy  method  which  is  similar  to  (6),  (7),  (8),  but  it 
requires  the  estimates  valid  in  the  half  space  and  in  the  exterior  domain.  He  note  that 
the  solution  decays  as  in  (1.10),  but  the  decay  rate  is  not  known,  cf.  (6|,  [7| ,  (8J.  in 


-3- 


.th*  following  wo  only  mention  the  argument!  for  the  Dirichlet  boundary  condition  (1.2), 
becauaa  the  Neumann  boundary  condition  (1.3)  for  6  can  be  treated  similarly. 


§2.  STATIONARY  SOLUTION 


Let  us  write  the  equations  end  conditions  for  the  stationery  solution  (p,u,  8): 


(2.1)  <puJ>  -  0  , 

J 

(2.2)  »  ♦  p  +  p*  -  (J(?  +  )  +  u’u*  611)  -  0  ,  i  -  1,2,3  , 

j  X1  Xi  j  i  \  j 


smsv 

KV?9x 


IMV  -4 
%>«u-  ~ 


( <9  )  -  ?  -  0 


(2.4) 


where  p  *  p(p,6)  etc.  The  stationary  problea  (2.1)-(2.4)  has  a  unique  solution  as 
La— a  2.1.  under  the  assumption*  A.1  ~  a. 3  there  exist  positive  constant  e  end  C  such 
that  if  ^  <  C,  i  ■  3  or  4,  the  problea  (2.1)-(2.4)  has  a  unique  solution  (p(x),0,'6) 
in  a  a— 11  neighborhood  of  (p,0,3)  in  (R2)3  satisfying 

(2.5)  Ip  -  plt  <  Cltlj,  4*3  or  4  respectively  , 
where  p(x)  is  determined  by  (1.8),  i.e., 

p(x)  p  (n,?) 

(2.6)  /  -  dn  +  *(x>  *  0 

I) 

p 

Proof.  Since  we  consider  a  eaell  neighborhood  of  ("p,0  ,"§)  in  (H2)3,  by  Sobolev’s  lemma 
we  may  suppose  |p  -  p|,  |u|,  |8  -  ?|  <  -j  min(p,"S).  Then  we  can  estimate  the  equalities: 

PpJn.'S) 

J  12.1)  K  J  -  dn  dx  -  0  , 

p 

(2.7)  /  (2.2)1  x  C^dx  -  0  , 

/  12.3)  *  (5  -  "Sldx  -  0  , 

where  (2.1),  12.2] A  and  (2.3)  danote  the  terms  on  the  left  hand  side  of  (2.1),  (2.2)1  and 


(2.3)  respectively.  Take  the  sua  of  (2.7)  and  integrate  it  by  parts  using  the  seen  value 
theorea  and  !«■■»  4. f.  Me  obtain  the  inequalityi 

idui2  +  id8i2  <  Cdopi  ♦  iui  ♦  |  e  -  'Si  .  ♦ 

1  c° 

♦  10  -  0I}(  I  DU  I2  +  ID0I2)  . 

« 

Therefore  if  lDpl,  lul^,  10  -  01^  is  small,  we  can  conclude 

(2.0)  u  -  0,  8  -  S  . 

If  we  substitute  (2.8)  in  (2.2),  we  have 


Ppt  n,  6) 

^  n 


p 


dn  +  #} 

Jx 

i 


0 


which  implies  (2.6) 


3.  LOCAL  AND  GLOBAL  EXISTENCE 

Let  us  rewrite  the  problem  (1.1),  (1.2)  by  the  change  of  variables 
(p,u,8)  •— *•  (p  +  p,u,8  +  9)  using  (2.6)  as  follows: 

uJP  +  (P  +  P)uJ  +  p  uJ  ■  0  , 


(3.1) 


(3.2) 


pt  +  ' x 


Xj  Xj 


u*  *  U^UA - - —  tutu1  +  u^  )  •*■  y'u*  + 

j  P  +  P  j  i  Tt  Xj 


Pp  Pg 

+  ~px  + 


pp  (p  +  p,e  +  9) 

(-P._  -  -  i )♦.  , 


P  +  P  Xi  p+pXi  (p  +  p)p  (P,9)  "i 


1.  <8  +  ?>p8  i  , 

*fr  Q  +  ■  as  —  — 

♦ 

v 


•  {( ic8  )  +  f}  . 

j  (P  +  P)c7  Xi  (P+P)cv  xj  xj 


(u,  9)  -  (u,  8)  | 

}Q 


(3.3) 


(P,u,8)(0,x)  »  (PQ,u0,80)(x) 


Further  we  rewrite  the  problem  (3.1 )— (3.3)  as  follows: 


(3.4)* 


(3.4)1 


(3.4)“ 


0  1-10 
L  (P,u)  5  p  +  uJP  +  -  f  , 

"  ^  XJ  Xj 


L1(p,u,8)  s  u*  -  yu*  -  (p  +  m')u3  _  +  P.P.  +  P,8  -  fi, 

c  *i*j  1  *L  * 


L4(u,9)  2  8—  <9  +  p  u^  -  f4  . 

^  J  ''j 


1  -  1 #2,3  , 


(3.5) 


(u,8)  I  ■  (u,  B)j  ■  0  , 

ho  I- 

(  P,u,  8)  (0)  -  (  P0fUQ,  0g)  , 


where  we  denote  the  constant  for  the  function  g  of  p  and  9  by  g  •  g(p,9),  and  also 
U  -  W/P.  M'  *  y'/P,  <  -  «/pc7,  p,  »  Pp/P,  P2  -  p^h  and  P3  «  "Spg/pc y.  The  terms  on 


■7' 


where 


*. 


Here  and  in  what  follows  we  do  not  write  ft  in  H*(ft)  and  L^CQ). 

We  will  obtain  the  global  solution  by  a  combination  of  a  local  existence  theorem  and 
some  a  priori  estimates  for  the  solution  in  X,  namely  that  for  the  norm  N . 

Proposition  3.1.  (local  existence) 

Suppose  the  problem  (3.1)-(3.3)  has  a  unique  solution  ( P,u, 8)  e  X(0,hiE0>  for  some 

h  >  0  and  consider  the  problem  (3.1)- (3. 3)  for  t  >  h.  Then  there  exist  positive 

constants  T , eQ  and  cQ ( cQ  /TT^<eo)  independent  of  h  such  that  if  N(h,h), 

1*1  <  e  ,  the  problem  has  a  unique  solution 

4  0 

(P,u,  8)  e  X(h,h  ♦  t  I  CflN(h,h) )  . 

The  proof  is  the  same  as  that  for  the  interior  problem  in  [8]  and  is  omitted.  Although  the 
local  existence  theorem  by  Tani  (10]  is  more  general,  we  need  it  in  the  form  of  Proposition 
3.1  to  extend  the  solution  globally  in  time  by  use  of  L2  energy  method. 

Propc sition  3.2.  (a  priori  estimates) 

Suppose  the  problem  (3.1 )— (3.3)  has  a  solution  (p,u,8)  e  x(0,h»E0>  for  given 
h  >  0.  Then  there  exist  positive  constants  and  C1  ( <  eQ,  <  EQ)  which  are 

independent  of  h  such  that  if  N(0,h),  1*1^  <  it  holds 

N(0,h)  <  CtN(0,0)  . 

If  Proposition  3.1  and  3.2  are  known,  the  global  existence  of  unique  solution  can  be 
proved  as  follows:  Choose  the  initial  data  ( , Uq  0Q )  and  the  potential  function  *  so 
small  that  it  holds 


Then  Proposition  3.1  with  h  “  0  gives  a  local  solution  (p,u,8)  e  X(0, TjCgNtO ,0 ) ) .  Since 
CgN(0 ,0 )  <  e1  <  Proposition  3.2  with  h  «  T  implies  N{0,  r)  <  C^tKO.O).  Then 
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Proposition  3.1  with  h  «  t  implies  the  existence  of  solution 


( p,u,8)  e  x<t,2tjc0n(t,t)),  e  x(o,2t»/ i  +  c*  n(o,t>)  . 

Hence,  since  /  1  +  Cg  N(0,T)  <  C 1  +  Cg  N(0,Q)  <  el ,  Proposition  3.2  with  h  « 
gives  N(0,2T)  <  C^N(0,0),  and  Proposition  3.1  with  h  “  2T  gives 

(P,U,8)  e  X(2T,3T,CgN(2T,2T)),  e  X(0,3t»/ 1  +  Cg  N(0,2Tj)  . 

Repetition  of  this  process  yields 
Proposition  3.3.  (global  existence) 

There  exist  positive  constants  £  and  C( ec  <  E. )  such  that  if  N(0,0),  141. 

U  4 

then  the  initial  boundary  value  problem  (3.1)-(3.3)  has  a  unique  solution 

(p,u,8)  e  x(0, •)CN(0,0 ) )  . 
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$4.  A  PRIORI  ESTIMATES 

First  we  recall  some  inequalities  of  Sobolev  type. 
Lemma  4.1. 


(4.1) 


(i)  Let  Q  be  any  domain  with  smooth  boundary.  It  holds 

0  <  o  <  1/2  , 

2  <  p  <  6  . 

(ii)  Let  (1  be  the  whole  space  R3,  the  half  space  R3,  or  the  exterior  domain  of 
a  bounded  region  with  smooth  boundary.  It  holds 


,fl  o 

<  CJfl 

C°(fl) 

H  (0) 

*f,L  (ft) 
P 

<  Clfl 

H  (0) 

(4.2) 


,f,L6(0)  <C,Df,L2(fl>  ' 


,f,L2(Q')  <C,D£,L2((1)  ' 


where  O’  is  any  bounded  subregion  of  (1. 

Proof,  see  for  example  [3],  [4]. 

Next  we  note  some  estimates  of  elliptic  system  of  equations  for  our  domain,  when  we 
regard  equation  (3.4)1,  i  =  1 .... ,4,  as  elliptic  with  respect  to  x  variables,  i.e., 

«  1  •  «  .  1  , 

WU  +  (M  +  M*  )uj  »  u  ♦  p.  p  +  p,  6  -f,  1  »  1,2,3  , 

xJxj  Vj  1  xi  P2 


(4.3)  x8  -  8  +  p  u3  -  f4  , 

xjxj  fc  3  Xj 

u|  «  u|  «  8 1  ■  81  «  0  . 

1 30  I-  I  an  I- 

Lemma  4.2.  Let  ft  be  the  half  space  or  any  exterior  domain.  He  have  for  k  »  2,3 

(4.4)  iD^ul  <  C{lutlk_2  +  lD(p,8)lfc_2  +  If  I  ♦  »ul)  , 

(4.5)  10*81  <  C{l8t<k_2  +  *D'i'k_2  +  ,f,k_2  +  ,D0,)  • 

The  first  estimate  is  well  known,  e.g.  (1).  The  last  Lz  norm  is  contained  on  the  right 
hand  side  because  of  the  unboundedness  of  our  domain.  The  second  one  is  given  in  [5] .  It 
contains  the  Lj  norm  of  the  first  derivative  as  the  last  term  on  the  right  hand  side  and 
is  stronger  than  that  containing  the  Lj  norm  of  the  function  itself  such  as  (4.4). 
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The  last  estimate  for  an  elliptic  system  concerns  Stokes  equation  in  8  which  comes 


from  (3.4)i,  i  -  0,...,4. 

Pu^  -  h  , 

i 

(4.6)  -Uu*  +  p  p  -  gL,  i  -  1,2,3  , 

J  j  i 

u|  -  a,  u|  -  0  . 

I  an  I- 

T  g*  4.3.  For  k  -  2,3,4  it  holds 


(4.7)  IdNiI2  +  lDk_1pl2  <  C{lhl2  +  Igl2  +  lal2  +  iDul2}  , 

2  Hk_V2<38) 


where  the  last  ten  on  hhe  right  hand  side  is  necessary  in  the  case  of  exterior  domain. 

Proof,  in  the  case  of  half  space  0  -  lemma  is  proved  by  Solonnikov  (9)  and  Cattabriga 

k—  1  /2 

[2],  where  the  definition  of  the  apace  H  (38)  is  also  given,  in  the  case  of  exterior 
domains  Finn  [3]  and  Heywod  [4]  obtain  lemma  provided  h  -  0.  Thus  we  only  need  a  slight 
improvement  for  h  4  0.  Let  us  extend  the  function  h  6  h’iSI)  to  a  function  h  e  H1  (*3) 
with  the  estimate 


Ihl  ,  <  Clhl 

H  (*  )  H  (8) 

Then  we  approximate  h  by  H1  functions  hn  with  compact  support,  i.e., 

h  ♦  h  in  h'i*3),  as  n  ♦  ■  . 
n 

Thus  we  can  define  the  function 


1  ,  hn<y> 

v*’  ■  j  /,  irrrr  * 


It  is  well  known  by  use  of  Fourier  transform 


ID2*!  .  ,  <  C»H  I  ,  , 

"  H^*3)  "  H^*3) 


Thus  we  have 
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ID'S  I  .  ,  *D4  L  ,  «•  V  «  Clh  I  , 

"h1!!))  “  V°  *  "hN*3) 

where  8*  is  any  bounded  subdomain  of  8  (cf.  (4.2)).  Since  C  is  independent  of  n, 
as  the  limit  of  n  ♦  •,  we  have 


*0  41  .  ,  ID4L  . qi  .  <  Clhl 

h’(8)  L2lQ  1  h\Q) 


u  -  V4  -  v  and  q  -  p  +  h/p1 P 


Now  we  can  put 


which  satisfy 

div  v  »  0 
*  / 

-U *  P,Vq  -  q 

and  reducss  to  the  above  case,  in  fact  we  have 

,DW,L2(0)  <  "*'^(8)  +  ,D^'l2(Q)  «  to“,L2<Q,  *  C,h'H1(n) 

and 


3/2  -  ■  ’+  •  <CID4I  <Clhl  , 

1 38  H'OO)  38hJ/Z(38)  HZ  ( 8* )  B 


(0) 


where  the  boundary  of  Q*  contains  38. 

Now  we  begin  to  obtain  the  energy  estimate  for  solution  of  equation  (3.4)1, 
i  -  0, . . . ,4,  with  (3.5). 

La— ui  4.4.  He  have  for  1-0  and  1 

I3*(p,u.9)(t)l2  +  f  lD3*(u,9)(s)l2  +  I32  TT  (s)  l2da 
t  0  t  t  dt 

(4.8) 

<  C{I32(  P,u,  8)  (0 )  I2  ♦  j  |A,|  «■  iaVl2ds> 

-  #1  •  t 


where 
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A  »  /  —  p(f°  -  u^P  )  +  uifi  +  ~  8f4  dx  , 
0  P  3  P3 


S  -  /  r1  Vf°  "  uJpxj)t  *  utft  +  5^  etftdx  * 


dp  .  j  .o-J 

TT  =  P«.  ♦  uJp  -  f  -  pu  . 
dt  t  x  x 


Ha  hava  also  for  k  «  0  and  1 


(4.10) 


lD3*(u,8)(t)l2  +  /  13*  (P,u,8)(s)l2ds 

c  0 


<  c{lD3*(u, 8)  (0 )  I2  +  I3*p(0 )  I2  ♦  I3*p(t)l2 


+  J  1d3*(u.8)(s)  I2  +  I3*(f°  -  ujp  )l2  ♦  I3*f(s)l2ds) 

0  *  j 


Proof.  Compute  tha  Integral 


/  /  —  pu.0  -  f°)  +  tt1a1  -  f1)  ♦  8(l4  -  f4)dxdt  -  o  . 

Oflp  P3 

Integration  by  parta  using  tha  boundary  condition  gives 

■jJ^p2+  |u|2+  —  82dx  +/  J  w|Du|2  ♦  (|i  ♦  p’ >(uj|  )2  ♦  —  x|  D  8|  2dxdt 


Xj  P3 


I  /  ?  P0  +  |u0|2  *  ?  4  l  V*  ' 

p  Vi  0 


where  Aq  Is  defined  by  (4.9).  If  we  use  the  notation  dp/dt  In  (4.9),  we  can  obtain 
(4.8),  t  •  0  from  this  equality.  The  time  derivative  can  be  treated  similarly,  because 
it  has  the  same  boundary  conditions.  Next  compute  the  integral 

/  /  p.  <L°  -  f°)  ♦  u^L1  -  fS  +  8  (L4  -  f 4) dxdt  •  0 
0  0  C  C 

Integration  by  parts  gives  by  use  of  Schwarz  inequality 


r 

* 


/  W|Du|2  +  (u  +  W*)(u^  )2  +  *e|  D0| 2  +  p  pu^  dx  +  1  /  /  p2  +  |u  |2  +  0?<5xdt 

fl  xj  1Xj  200t 

<  /  U | DUg  | 2  +  (W  +  w',<uo,x  )2  *  le'D9o'2  +  p1p0U0,x  ^ 

t  i  -  i 

+  C  /  /  |D(u,8)(s)|  ♦  |f|  dxdt  , 

o  a 

where  f  “  (f°  -  u^p  .f1 ,f2,f3,f4).  If  we  use  Schwarz  inequality  for  the  tern  pu^  ,  we 

Xj  Xj 

obtain  (4.10),  i  ■  0.  The  estimate  (4.10),  1  ■  1  is  obtained  similarly. 

In  the  following  we  first  treat  the  case  of  half  space  Q  -  m2  which  is  simpler  to 

explain  the  way  of  obtaining  estimates  than  that  for  exterior  domains  treated  later. 

Since  the  tangential  derivatives  of  the  solution  of  (3.4)  satisfy  the  same  boundary 

conditions  (3.5),  we  can  obtain  the  estimates  for  these  similarly  to  the  above  lemma  4.4. 

Let  us  denote  the  tangential  derivatives  by  )■  ()  ,3  )  and  integrate  the  equality  for 

X1  x2 

each  k  -  1.2,3  by  use  of  integration  by  parts 

P1  »k,T0  _0.,k  4k,,i  ,1,Jt  l  _  . 

—  3(L  -  f  lip*  3  (L  -  f)3u  -  0. 

P 


Thus  we  have 

Lemma  4.5.  For  k  *  1,2,3 


(4.11) 


I3k(p,u)(t)l2  +  J  lD3ku(s) I2  +  I3k  (s)l2ds 
0 


<  C{l3k(p,u)  (0)  I2  +  /  I3kf°(s)l2  +  I3k_1f(s)l2  +  |v  ..Ids)  , 
0  K  1 


where 

(4.12) 


“  /  — ■  3kp3k(f°  -  u^Px  )dx  for  each  k  »  1,2,3 


0  p 


Then  we  have  to  obtain  the  estimates  for  the  normal  derivatives  of  solution.  To  do 
that  we  use  the  following  equations  from  (3.4). 
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(4.13) 


where  we  note  the  second  derivatives  of  u  at  the  last  two  terns  on  the  right  hand  side 

contain  one  tangential  derivative.  Multiply  (4.14)  by  p  and  (d(Vdt)  respectively 

X3  X3 

and  integrate  them  respectively.  He  obtain  after  integration  by  parts 


I  3M-+-iL  p*  dx  +  /  /  p1P^  dxdt 


2P  3  0 


/  t  P»  dx  +  /  /  2t>-j  ^  ■  (-u'l  p  +  2u3  p  )p  dxdt 

2p  u'x3  0  2p  *3  3  j  X3 


and 


+  /  /  {-u3  -  p  8  +  p(u3  +  u3  +  u1  +  u2  ) )p  + 

'  *  t  *2  X  XX  XX  XX  XX  X 

0  3  1  1  2  2  1  3  2  3  3 


♦  *  JL’  t®  +  f3)Px  dxdt 

P  3  3 


p.  t  t 

<  ~  {  !  P‘  dxdt  +  c{lDp(0)l  *  [  |B.  n  |dt  + 

2  0  X3  0 


+  /  In  I2  +  iDJul2  +  ID8I2  *  lDf°l2  +  If3  l2ds  , 
0  * 
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V 


♦  /  lu  I2  +  Id8I2  ♦  Id9uI2  +  lDf°l2  ♦  If2 l2dxdt } 

A  ^ 


respectively.  Thus  wa  have  obtained  the  following 
Lew  a  4.6.  For  k  +  i  ■  0,1,2  it  holda 


where 

(4.16) 


I3k9*+1p(t)l2  ♦  /  I3k9*+1  p(a)  I2  +  *^»*+1(|f)(a)l2da 
<  c{iDp<0)i^t  ♦  /*  la^’ajbui2  ♦  i3ka*uti2 
+  «0C«.WC»<rt*  ♦  «<»♦  lBk>l|da}  . 

*k«t  *  / 1***,4©.  ■  ' 


and  here  the  sussution  ie  not  taken  for  k  and  i. 

Proof.  For  k  ♦  t  •  1  and  2  we  differentiate  the  equation  (4.14)  by  9k92  and  aultiply 

it  by  9k9*+1p  and  3k3*+1  respectively,  integration  by  parts  gives  (4.15)  and  (4.16) 

in  this  case  similarly  to  that  for  k  ■  t  -  0. 

Last  we  use  Issbm  4.3  for  Stokes  equation  (4.6)  with  ul  -  0,  where  h  and  gl 

( 90 

have  the  following  explicit  forms. 
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(4.17) 


h  -  f°  -  & 
dt 


1  “  *  -  J"*1'  hx  -  p20x  +  f1.  i  -  1,2.3  . 


P  i 

4.7.  For  k  ♦  i  -  0,1,2  we  have 


lD2+t9kul  ♦  lD1+t9kpl 


(4.18) 


<  c{.utVt  «■  I9k(|f )lm  -  «08.  ♦  'f°'1+k+t  ♦  '*'*♦*}  • 


Now  we  can  combine  the  above  lemmas  4. 1-4.7  to  obtain  necessary  a  priori  estimates. 
First  we  obtain  the  H  version  of  norm  N(0,t),  i.e., 

•  p,u,  8(t)  <2  +  lpt(t)l2  +  'vVt)|2 


(4.19) 


+  /  lpt(s),Dp(s)l2  +  lu  ,8t(s)l2  +  «D(u,6)(s)l2  +  (e)l2ds 

<  c{lp,u,8(0)l2  +  sup  0f°  -  u3p  l2,lf(s)l2} 
o<s<t  3 

+  f  if°i2  +  «f°i2  +  ifi2  ♦  iuJp  i2  +  i  i\ i  ♦  i  \\  .ids . 

n  a  *  ■  **  *  l_n  *  k+t*0 


It  is  proved  by  fifteen  steps  as  follows: 
(i)  By  lenma  4.4,  1  -  0 ,  we  have 


•p,u,8(t)l2  +  /  lD(u,8)(s)l2  +  t~£  (s)l2ds 

0 


(4.20) 


<  c{lp,u,8(0)l2  +  /  lf°(s)l2+  | A  I ds } 

0 


(ii)  By  lemma  4.4,  k  “  0,  and  (4.20)  we  have 


lD(u,9)(t)l2  +  /  lp.,u.,8  (e)l2ds 

A  fc  &  V 


(4.21) 


< . +  C  { lD(u,  8)  (0 )  I2  «■  lp(0)l2  ♦  /  lf(s)l2ds)  , 

0 
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where 


means  the  terms  already  appeared  on  the  right  hand  side  of  Inequalities  (4.20)- 


(4.32),  especially  (4.20)  for  (4.21). 

(Ill)  By  lenma  4.4,  1  «  1,  we  have 


,pt*vVt,,a  +  f  "’'vV"1'2  +  ,^t<B),2ds 


(4.22) 


<  c{«Pt,ut,9t(0)l2  +  /  lf®(s)l2  +  1*^  Ids } 


(iv)  By  leans  4.2  for  0  and  by  (4.20)-(4.22)  we  have 


(4.23)  lD20(t)l2  +  /  ID20(s)  l2ds  < 

0 


+  C{lf4(t)l2  ♦  /  lf4(s)l2d*} 

0 


(v)  He  have  the  estimate  of  tangential  derivatives  of  p  by  leans  4.S,  k  ■  1 


•3p(t)l2  +  /  lo3u(s)»2  ♦  I0(^)(s>»2ds 


(4.24) 


« . +  I3p(0)  I2  ♦  C  /  13f°(s)l2  +  lf(s)l2  +  |A  Ids  . 

0 


(vi)  The  estimate  of  normal  derivative  for  p  and  — §  follows  from 

at 

k  +  t  -  0  and  from  (4.24),  (4.21)  and  (4.20) 


lemma  4.6 


(4.25) 


I3,p(t) I2  ♦  /  I3_p(s)l2  +  13,  (s)l2ds 

3  o  3  3  dt 

< . +  c{l3  p(0)l2  ♦  /  lf°l2  Ifl2  +  |B  Ids} 

J  1  UpU 


(vii)  Than  we  have  the  second  derivative  of  u  by  lenma  4*2  for  u  end  k 
by  (4.19),  (4.21),  (4.22),  (4.24),  (4.25). 

(4.26)  l02u(t) I2  <  .  ♦  Clf(t)l2  . 


2,  and 
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(viti)  Further  since  e  L.^tiH1),  we  heve  by  lean*  4.7,  k  ♦  i  «  0  end  by 
at  2 

(4.21),  (4.20),  (4.24)  end  (4.25) 


(4.27)  /  lD2u(»)l2  +  lDP(e)  l2de  < . +  /  If0  I2  +  lfl2de 

0  0 


(ix)  By  lean*  4.2  for  6,  k  »  3  end  by  (4.27),  (4.21)  end  (4.22)  we  heve 


(4.28)  /  lD39(«)l2da  < . +  /  lf4(s)l2ds  . 

0  0 


(x)  By  leans  4.5,  k  «  2  we  heve 


(4.29) 


(xi) 


(4.30) 


I32p(t)l2  +  /  ID22u(e)  I2  +  I32  (a)  l2da 

0  * 

<  c{ia2(p,u)(0)l2  ♦  /  ia2f°(e)l2  ♦  I3f(e)l2  ♦  |A  |de}  . 
0 

By  Inane  4.6,  k  -  1,  i  -  0  end  by  (4.29),  (4.22)  we  heve 

iaa3p(t)i2  ♦  Jfc  iaa3P(s)i2  +  iaa3  f*  (a)i2d* 


< . ♦  C{l33,p(0)l2  +  /  lf°(s)l2  If  1?  ♦  |B.  .Id*}  . 

3  0  *  1  1,0 

By  Leans  4.7,  k  «  1,  l  -  0  end  by  (4.29),  (4.30)  end  (4.22)  we  heve 

(4.31)  /  lD2au(*)l2  ♦  lD3p(s)l2da  < . +  /  lf°(e)l2  ♦  If ( s)  l2da  . 

0  0 

(xiii)  By  leeaw  4.6,  k  »  0,  i  “  1  end  by  (4.31),  (4.22)  we  have 

(4.32)  I32p(t)l2  +  /  I32p(a)l2  I32  (a)l2d*  < . ♦  c{lp(0)l2  +  f  |B.  ,|da} 

J  0J  **  at  2  Q  0  *  1 
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(xiv)  By  It 


k 


4.7,  k  -  0,  t  -  1  and  by  <4.29),  (4.30),  (4.32)  and  (4.22)  we  have 

(4.33)  J*  I03u(a)l2  +  lD2p(a)l2da  < . 

0 

(xv)  Then  by  the  equation  (3.4) 

(4.34)  lp„(t)l2  ♦  /  Ip  (a)  l2da  < . +  If0  -  uV  I?  ♦  /*  If0  -  uJp  (a)  if  da  . 

t  i  o  *  1  xj  1  o  j  1 

Thua  if  «e  aua  (4. 20)- (4. 34),  we  arrive  at  the  H2  energy  eatiaate  (4.19). 

To  elevate  the  differentiability  once  to  obtain  the  eatiaate  of  non  N(0,t)  we  can 
repeat  the  above  argument  beginning  froa  lane  4.4,  k  -  1  and  by  one  of  leaaw  4.2, 
k  -  3,  lama  4.5,  k  -  3,  lean  4.6.  k  +  t  -  2  and  lean  4.7,  k  +  i  -  2.  Therefore 
we  arrive  at  the  eatiaate  for  N(0,t). 


11(0, t)2  5  lp,u,9(t)  I2  ♦  IPt(t)l2  ♦  (h»t,et(t)i2 


♦  J  IPt.DP(a)l2  +  lut,9t(a)l|  ♦  lD(u,  6)  (a)  l2da 


(4.35)  <  c{lp,u,9(0)l2  +  aup  {If0  -  u3Px  (a)  I2,  If  (a)  I2 


0<a<t 


♦  /  l«°(a)l2  *  If  (a)  I2  ♦  If  (a)  I2  ♦  l(ujp.  ).  I2  +  hijp  I2 

0  ■*  '  1  *j  ‘  *  j  2 

4  2 

♦  I  lA-l  ♦  [  1 8.  . | da }  . 

k-0  *  k+twO  *'* 


Laat  wa  have  to  ahow 
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i 


sup  {f°  -  u3Px  (s)l2,lf(s)l2>  +  /  lf°ts)l*  + 


0<s<t 


(4.36) 


+  Kf°  -  u3 P  ).<s)l2  +  lf.(s)l2  +  lf(s)|2  +  I  |A.|  *  l  |B.  ,|ds 

xj  '  c  *  jpmQ  K  «_n  *«* 


<  C2(U(0,t)  +  UI4)N  (0,t) 


It  is  proved  by  use  of  leans  4.1  end  integration  by  parts.  We  show  only  the  ten  Aq  and 
oeLit  the  proof  of  the  other  terms  which  can  be  treated  similarly.  Let  us  recall  (4.9)  and 
compute  the  following 


I  /  P(f°  -  u3p  )dx|  -  |/  p{<  p  -  p)u3  -  pu3}  dx|  -  |  /  p  {( P  -  p)u3  -  pu3 )dx| 

j  j  j 

<  iDpI  {(/  (p  -  p)2|u|2dx]1,/2  +  (/  p2|u|2dx)1//2 } 


<  IDpI { Ip  -  pi  lul  +  Ipl  lul  } 
3  L6  “3  6 


<  Cl DP I iDul ( I • «3  +  Ipl,)  <  C  N(0,t)2{l*l3  ♦  N(0,t)> 


1/  ui(( — -  u)ux  +  — — ;  ux  ux  }dx| 
P+P  )j  P  +  P  )j 


<  /  \\  ( — ^~z  -  >0ux  I  +  • ( — 3 — z)x  ux  |d* 
)  P  +  P  J  P  ♦  P  j  J 


<  | - K—  _  y)  lDu,2  +  l(— I  lul.  IDUI 

P+P  °  P  +  P  *3  L3  L6 


<  C  N(0,t)  (Ulj  +  N(0 , t)  ) . 


The  remaining  terms  in  aq  can  be  treated  in  the  same  way  as  above 


Finally  let  us  turn  to  the  case  of  exterior  domains.  In  this  case,  since  we  cannot 
generally  designate  a  coordinate  system  all  over  Q  the  directions  of  which  are  consistent 
with  the  normal  and  tangential  directions  on  the  boundary  312  as  in  the  case  of  half 
apace,  we  have  to  modify  the  lemmas  4.5  through  4.7.  To  do  that  we  shall  separate  the 
estimates  of  the  solution  into  that  over  the  region  away  from  the  boundary  2(1  and  that 
near  the  boundary  3Q.  Let  Xg(x)  be  any  fixed  cut-off  function  in  C  ( ft)  such  that 
support  Xg  C  (1  and  Xg  =  1  outside  of  a  bounded  region  O' .  Then  we  have  the  following 
as  the  estimate  on  the  region  away  from  the  boundary. 

Lemma  4.9.  For  k  »  1,2,3,  it  holds  that 


•XgD  P<t)r  +  /  lx0DKp(s)rds 


(4.37) 


<  cjlDkp(0)«2  +  lDk-1u(0)l2  +  IDk~^u(t)  l2 


+ 1  ,f°  -  +  '<i  +  ,Dtu'9),k-i  +  ' 

IX0Dlt(P,u)(t)i2  +  t  »x/+1u(s>.2  +  ly>k  (s)l2ds 


(4.38) 


<  c{lDk(p,u)(0)|2  +  /  lDk_1fl2  +  lDkf°l2 
0 


where 


♦  «X0DkP«2  *  ID(u,8)I2_1  ♦  jA^lds} 


0  P 

and  here  the  summation  is  not  taken  for  k. 
Proof,  compute  the  integral 
t 


\+1  -  /  ~  Xg (Dkp)Dk( f°  -  ujPx  )dx  , 


/  /  liLl-Hi  .  f0,  p  +  ,2(Li  .  fi)p 

0  Q  p  i  i 
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Altar  integration  by  parts,  we  have 


/  xjlDPl2  ♦  )UU1PX  dx  |  +  /  /  p1  )L|Dp|2dx<U 

fi  2P  i  '0  0  Q 

•  l  l  -a«b*b.*j,ijp*i  +  2^Xo,Sci“xjexj  +  4‘\ 

-  *i4\\  *  (2X0*0.xiui  +  X0Uxj)(u3pxJ  -  f0  +  pu2j)dx  +  S 


da 


Since  the  support  of  DXg  is  contained  in  ft' ,  we  can  sake  use  of  (4.2)  and  consequently 
Schwarz  inequality  give*  (4.37),  k  »  1.  The  other  cases  k  *»  2,3,  are  obtained 
eiailarly.  To  prove  (4.36),  coapute  the  integral 


/  /  xf  —  (00)D(L° 

0  Q  p 


f°)  ♦  Xq(DU^)D(I*^  -  f1)dxds  -  0 


Integration  by  parts  gives 


X„  r-  |DP|2  ♦  XplOu|2dx  )fl  +  /  /  j£»*|D2u|2 

’  I  l  &2(Dux1,D8  *  2\>\>,x/2«Du1>09 

+  2>fe)t0.xiP1<Dtti,Dp-  X°VuKjiej  '  2%\),xj 


+  )L(U  +  MMIDu^ 

1 


fV  dx  +  A^ds  , 


|2dxds 


which  ieplies  (4.38),  k  ■  1  after  using  Schwarz  inequality.  The  cases  k  •  2,3  are 
obtained  similarly . 

Next  let  us  establish  the  estiaates  near  the  boundary.  To  do  that  we  choose  a  finite 
nuaber  of  bounded  open  sets  {0^}^  in  X  such  that 

M 

U  0,  3  »0  , 

3-1  3 

and  in  each  set  0 ^  we  choose  local  coordinates  (f, +,r)  as  follows; 


(i)  the  boundary  0^  n  3n  is  the  image  of  smooth  functions  yA  “  v 1  (  4,  4)  satisfying 
(ax.,  taka  the  local  geodesic  polar  coordinate) 

(4.39)  ly^l  -  1,  yjyj  -  0.  ly^l  >  4  >  0  . 
where  4  is  some  positive  constant  independent  of  J  •  1,...,N. 

(ii)  any  x  in  0 j  is  represented  by 

(4.40)  X*  -  xi(<i,4,r)  -  rni(<i,4)  ♦  yi(*,4)  , 

where  nN  ♦» ♦)  represents  the  external  unit  normal  vector  at  the  point  of  the  boundary 
coordinated  (9,f).  Here  and  in  what  follows  we  omit  the  suffix  j  for  simplicity.  Let 
us  define  the  unit  vectors  e*  and  e*  by  e*  “  y*  and  e*  ■  y^ly^l.  Then  Prenet- 
Ser rat’s  formula  gives  that  there  exist  ssooth  functions  ( a,  B,  Y,  a* ,  S' ,  y’ )  of  (9,4) 
satisfying 


An  elementary  calculation  shows  that  the  Jacobian  J  of  the  transformetin  (4.40)  is  given 
by 

(4.41)  J  •  Ix^  «  x^|  -  |y^|  +  (afy^|  ♦  8')r  t  («  ('  .  8a')r2  . 

By  (4.41)  we  can  see  the  transformation  (4.40)  is  regular  choosing  r  small  if  needed. 


Therefore  the  functions  (9,#,r)  (x)  make  sense  and  is  calculated  as 

i 

\  ’  1  <  V  Vi  “  ?  +  • 


(4.42) 


\  "  i  <xr  *  Vi  "  J  <C*}  *  "i'  ’ 


rx1  ■  J  ,x4  *  Vi  “  ni  - 


where  A  -  |y^|  +  8'r,  a  ■  -ro',  C  -  -Sr,  D  «  1  ♦  or  and  J  -  ad  -  BC  >  0.  Hence 
(4.42)  implies 
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«. 


8*a  “  J  (Aet  +  Be2,%  +  i  (C*t  *  “•a’3*  +  n^3r  * 


Thus  in  sach  0 j  we  can  rawrita  the  aquations  {(3.4)i}A_Q  in  tha  local  coordinates 
as  follows: 


L°  “  dt  +  J  <(A*i  +  B*2,ut  +  <Ce*  *  °«2)u8  +  Jnlurl  “  f°  • 

iZ-uf  ^II^bV  . -  1  ‘  -2  ‘  -2‘  1  ■  -2  i 


fc  J2 


w  ♦  2 (AC  +  BDju^  +  (C  +  D  )u^  ♦  J  Urr)  ♦ 


*  first  order  terns  of  u  and  8 

*  ;  <“!  *  *  p,»),  *  3  •  p,»),  * 

*  .  iii: ,« 

P  Pi 

•hare  wa  note  that  J2  -  (AC  +  BD)2  -  (a2  +  B2)(c2  +  D2).  Let  us  denote  the  tangential 
derivatives  by  3  -  (3^,3^)  as  before  and  let  JCj <  1  <  J  <  M)  be  any  fixed  cut-off 
function  in  C*( 0^ ) .  Estimating  the  integral  for  k  “  1,2,3 

/  /  X?  —  <3kp)3ka°  -  f°)  ♦  )^(8ltu1)3k(l.1  -  f1)  J  dfMdrds  -  0 

0  0  3  p  ^ 

in  the  similar  way  as  in  leonas  4.5  and  4.9,  we  can  get 

4.10.  For  any  positive  c  and  k  “  1,2,3,  it  holds  that 


IXj3k<P.u)(t)l2  +  /  Ix^uls)!2  +  lxj3k  (s)l2ds 


<  c{lxj3k<P,u)(0)l2  +  /  elDplj^,,  ♦  (1  +  e“1)lD(u,8)l2_1 


+  ,f0,£  +  ,f,:-i +  I'w*-)  * 


where  A  *  /  —  X?(3kp)3k(f°  -  u3p  )dx,  and  here  the  summation  is  not  taken  for  k. 
3'*M  OP3  Xj 

In  order  to  estimate  tha  normal  derivatives  we  make  use  of  the  equation 


3r<L0  -  f°)  ”  0  and  nSl.1  -  f*)  ”  0  which  have  tha  form 
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©r  +  J  {(Ae1  +  Be2)ur*  *  (Ce}  +  °*2)u^  +  ^"rr* 

0 

+  first  order  terms  of  u  *  , 

* 

(4.43)  "^“t  +  ^2  +  B^n^u^  +  2(AC  +  BD)u^  +  (C2  +  02)niU^  +  J2n^u^r) 

J 

+  first  order  tents  of  u  and  6 

^©r  +  Plpr-n 

p  p 

Eliminating  n^u^r  from  (4.43),  we  get 


2u 


^dt^r  *  p1Pr  “  “ni“t  +  ^2  *(a2  +  *  2<AC  +  BDlnNi*^  + 


(4.44)  +  (c2  ♦  D2)niu^  -  J(Ae*  ♦  Be^Ju^  -  J(Ce*  *  De^u^J  + 

i  i  u  ♦  u*  0 

+  first  order  terns  of  u  and  B  +  n  f  +  — r-  f^  . 

P 

If  we  apply  3k3*  (k  +  l  "  0,1,2)  to  (4.44),  multiply  it  by  ** )r  and 
and  integrate  them  respectively  in  the  similar  way  as  in  lemmas  4.6  and  4.9,  we  can  have 
Lemma  4.11.  For  k  +  t  «  0,1,2  it  holds  that 


where 


•x^a^W)!2  ♦  /*  iXj3k3*+1p(s>i2  +  »xjaka*+1(|^)(8)i2ds 

<  (c{lDp(0)l^.t  +  f  Ix^VduI2  +  lXj3k3*utl2  + 

+  .D(u,e)l2>j{  +  .f0.2^,  ♦  .f.J+t+  lBj>k/tlds>  , 

B  -  /  X2(3k3*+1[^)  "  3k**+1p  )3k3*+1pdx  , 

“j,k,l  „  3  r  ldt'  r  V  r  **“  ' 


and  here  the  summation  is  not  taken  for  k  and  t. 
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Last  we  have  to  get  a  leasts  corresponding  to  leima  4*7  •  Evidently  the  statenent  of 

leasts  4.7,  k  ♦  l  ■  0  holds  also  for  the  case  of  exterior  dooains,  if  we  add  the  ten 

iDul  on  the  right  hand  side  of  the  Inequality,  (cf.  lemma  4.3).  Kext  operating  3  , 

k  -  1,2  to  Stokea  equation  (4.6)  with  (4.17)  and  u|  “  0,  we  have 

30 


(4.45) 


p(y“ul>Xi  ■  yk<°  -  x3ak(ff)  ♦  pxj^^u1  . 

+  p1(xj3kp)*1  "  "  p2X33>texi  *  *1<J,«1*P  ‘  MX3.xtxta 

-  *  y*'1  -  ^  ykdf) +  ^  ' 


^3  u 


3  u  ”  0  . 

30 


Thus  we  can  apply  lemma  4.3  to  (4.45)  and  consequently  we  have 

Lemma  4.12.  For  i  “  0,1,2,  it  holds  that 

,_2+i  .  .  ,_1+i  .  .  „  ^  -dp,  ^  „0, 


ul  +  lD1+tpl  <  c{lut«t  +  ♦  ,f0,1  +  t  ♦  ,f,t  ♦  lD(u,e>l}  , 


and  for  k  «  1,2,  t  *  k  »  1,2,  it  holds  that 

-  ,.2+t-k  ■  „1+i5k  ■ 

IXjD  3  ul  +  Ix^D  3  pi 


<  c{,Vk+* +  *yk  S'i+* ♦  ,f0w  +  ,fw  ,dpWi  *  ,D(u'e,w 

Combining  the  above  lemmas  4. 1-4.4  and  4.9-4.12  in  the  similar  way  as  in  (i)-(xv)  for 
the  half  space,  we  can  get  the  exactly  same  inequality  (4.35)  and  the  same  lemma  4.8  if 


only  replacing 


l  l*J  +  I  |Bk  ,1 

k-0  k+t-0  ' 


N  4 


by  |A  I  +  |A  I  +  l  |A*|  +  [  {  l  I*  |  +  l  |B  1} 

k-2  J-1  k-2  3,K  k+twO 


Thus  the  proof  of  Theorem  1.1  is  completed. 


umncis 


II]  A;pon,  S.,  oouglis,  A.  and  Nirenberg,  L. ,  estimates  naar  tha  boundary  for  solutions  of 
alliptic  partial  diffarantial  aquations  satisfying  ganaral  boundary  conditions  II, 
Coam.  Pura  Appl.  Math.,  J7.  (1964)  3S-92. 

[2]  Cattabriga,  "Su  un  problaaa  al  contorno  ralativo  al  sistaaa  di  aquasioni  di  Stokes, 
Rand.  Mat.  Sem.  Univ.  Padova,  31,  (1961)  308-340. 

[3]  Finn,  R. ,  On  tha  axtarior  stationary  problaa  for  tha  Naviar-Stokas  equations,  and 
associatad  perturbation  problaaa,  Arch.  Rat.  Mach.  Anal.,  J9_  (1965)  363-406. 

[4]  Haywood,  j.  G. ,  A  uniqueness  theorsa  for  non-stationary  Havier-Stokes  flow  past  an 
obstacle,  Ann.  Scuola  Norm.  Sup.  Pisa  IV  6,  (1979)  427-445. 

[5]  Lady shanskaya ,  0.  A.,  The  Mathematical  Theory  of  Viscous  Incompressible  Flow,  Gordon 
and  Breach,  Haw  York,  1969. 

16]  Matauauira,  A.,  An  energy  method  for  tha  aquations  of  notion  of  compressible  viscous 
and  heat-conductive  fluids,  Univ.  of  Wlaconsin-Madison,  HRC  Technical  Sunnary  Report 
•2194,  1981. 

[7]  Matsunura,  A.  and  Hishida,  T.,  The  initial  value  problem  for  tha  equations  of  motion 
of  compressible  viscous  and  host-conductive  fluids,  Proc.  Japan  Acad.  ser.  A,  55 
(1979)  337-342. 

[8]  ......  Initial  boundary  value  problems  for  tha  aquations  of  motion  of  general  fluids, 

to  appear  in  Proc.  of  5th  Internat.  Symp.  on  Computing  Methods  in  Appl.  sci.  and 
Bngin.,  Dec.  1982,  IMRIA,  Versailles,  Franca. 

[9]  Solonnikov,  V.  A.,  A  priori  estimates  for  certain  boundary  value  problems,  Sov.  Math. 
Dokl. ,  2.  (1961)  723. 

[10]  Tani,  A.,  On  tha  first  initial-boundary  value  problem  of  compressible  viscous  fluid 
motion,  Publ.  RXM8,  Kyoto  Univ.,  U.  (1977)  193-253 


AM/TN/ed 


-29- 


SECURITY  CLASSIFICATION  OF  THIS  FACE  (Wkan  Dot,  Holered) 


|  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

I.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (end  Subtitle) 

INITIAL  BOUNDARY  VALUE  PROBLEMS  FOR  THE 

EQUATIONS  OF  MOTION  OF  COMPRESSIBLE  VISCOUS 

AND  HEAT-CONDUCTIVE  FLUIDS 

S.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Summary  Report  -  no  specific 
reporting  period 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AuTHORfcJ 

Akitaka  Matsumura  and  Takaaki  Nishida 

S.  CONTRACT  OR  GRANT  NUMBERfoJ 

DAAG29-80-C-00  41 

».  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Mathematics  Research  Center,  University  of 

610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 

Work  Unit  Number  1  - 
Applied  Analysis 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U.  S.  Army  Research  Office 

P.O.  Box  12211 

Research  Triangle  Park,  North  Carolina  27709 

12.  REPORT  DATE 

September  1982 

IS.  NUMBER  OF  PAGES 

29 

14.  MONITORING  AGENCY  NAME  4  AODRESSflf  dtllerent  from  Controlling  Otttca) 

IS.  SECURITY  CLASS,  (ot  tble  report) 

UNCLASSIFIED 

ISa.  DECLASSIFICATION/DOWNGRAO'NG 
SCHEDULE 

IS.  DISTRIBUTION  STATEMENT  (ol  Me  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  tho  abetted  entered  In  Block  20.  II  dlllerent  from  Report) 


IS.  SUPPLEMENTARY  NOTES 


IS.  KEY  WORDS  (Continue  on  roretee  aids  II  naesaaar y  and  Identity  by  Moc k  number) 

Equation  of  general  fluids.  Initial  boundary  value  problem, 
Global  solution  in  time 


20.  ABSTRACT  (Continue  on  rororoo  aide  II  neceeeery  and  Identity  by  block  number) 

The  equations  of  motion  of  compressible  viscous  and  heat-conductive 
fluids  are  investigated  for  initial  boundary  value  problems  on  the  half  space 
and  on  the  exterior  domain  of  any  bounded  region.  The  global  solution  in 
time  is  proved  to  exist  uniquely  and  approach  the  stationary  state  as  t  -*•  «*, 
provided  the  prescribed  initial  data  and  the  external  force  are  sufficiently 
small. 


00  ,  jSTn  1471  EDITION  OF  I  NOV  SS  IS  OBSOLETE  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAOE  (*Mr  85»  Smtorad) 


